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method of measuring and correcting the distortion in a 
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channel distortion. Lagrange interpolation of the 
received values of the known pilot symbols is used to 
estimate the channel distortion within the pilot symbols. 
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Description 

Held Of The Invention 

* [0001 ] This invention is related to digital signal processing, particularly to compensate for channel distortion in dig- 
ital mobile radio transmission. 

Background Of The Invention 

10 [0002] Digital mobile radio communication is plagued by distortion in the transmission channel. The Doppler effect 
due to a moving transmitter or receiver and multiple paths of propagation due to reflection are two principal causes of 
this cf stortion. Compensating for this distortion is particularly important in a Direct Sequence Code Division Multiple 
Access (DS-CDMA) system, Accurate channel distortion compensation enables the DS-COMA system to operate in a 
coherent delay-locked tracking loop (DLL) mode which has numerous advantages over a noncoherent delay-locked 

15 tracking mode. Coherent DLL mode of operation reduces the background noise by as much a 3 dB compared to the 
noncoherent DLL mode and, for the same BER (Bit Error Rate), requires less transmitted power output 
[0003] Distortion in a digital mobile radio channel is predominately manifested as amplitude and phase variation of 
the channel gain. If the channel gain can be accurately estimated, then good compensation for the channel distortion 
may be achieved. One approach to obtaining a good estimate of the channel gain is based on insertion of known synv 

20 bote at regular intervals in the data stream. When these symbols are recovered in the receiver any deviation from their 
known values is taken to be caused by the time-varying channel gain. From these deviations, the channel gain may be 
estimated. Various methods are known for such insertion and recovery of known symbols, but all are subject to limita- 
tions. For example, one such estimation method uses only a few symbols at a time which makes this method vulnerable 
to noise. Another exemplary method uses an experimental method to determine a fixed set of estimation formulae. 

25 These fixed estimation formulae do not adjust to the changing condrtions that exist during operation in a real channel. 

Summary Of The Invention 

[0004] The present invention provides an efficient method of estimating and compensating for the rapid amplitude 
30 and phase fluctuations in a digital mobile radio transmission system based on the insertion of known pilot symbols peri- 
odically into the data stream. Multiple consecutive symbols are inserted and the consecutive received values averaged 
to overcome random noise. These averaged values provide samples of the complex, random channel gain. These sam- 
ples are used for interpolation by the Lagrange polynomial method to yield the values of amplitude and phase distribu- 
tion undergone by each channel data symbol. 

35 

Brief Description Of Tfrfi BflMEfiS 
[0005] 

40 Figure 1 is a block diagram of a typical DS-CDMA receiver. 

Figure 2 is a diagram of quadrature phase shift keying (QPSK) modulation 

Figure 3 is a diagram of a data stream divided into blocks with pilot symbols inserted at the beginning of each block. 

45 

Figure 4 depicts a basic channel estimation and compensation processor. 

Figure 5 depicts a channel estimation and compensation processor using the weighted multi-slot averaging 
(WMSA) method. 

so 

Figure 6 depicts a channel estimation and compensation processor according to the method of the invention. 
Detailed Bescrjpjlpj] 

55 [0006] The invention is directed to an improved method for determining channel gain in a wireless communications 
system where data is transmitted in a digital mode. Application of the methodology of the invention is described here- 
after in terms of a preferred embodiment based on a DS-CDMA wireless system. It will, however, be apparent to those 
skilled in the art that the inventive methodology may be applied for a variety of digital wireless systems. 
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[0007] Figure 1 is a block diagram of a typical DS-CDMA receiver system of the present art. An antenna 100 
receives a spread-spectrum radio frequency (RF) signal A typical carrier frequency is 2 gigahertz (GHz) with a typical 
bandwidth of 5 megahertz (MHz). An RF receiver/demodulator 102 converts down from the carrier frequency The 
result is an analog baseband signal, modulated by an encoded data stream, typically at 32 kilobits per second (kbps), 
5 and re-modulated by a spreading signal, typically at 4.096 megachips per second (Mcps). ("Chip" is the standard term 
for one cycle time of a spreading signal.) A matched filter 104 removes the spreading signal by correlation with the 
phase of the spreading signal. The remaining analog baseband signal is processed by a channel estimation and com- 
pensation processor 106 to produce an estimate of the encoded data. 

[0008] Because the RF signal received at antenna 100 will generally include the superposition of multiple images 
to of the source transmission - representing different paths of propagation with different delays -- multiple instances of the 
filter 104 and compensation processor 106 may be used to extract an estimate of the encoded data stream from each 
of several of the strongest paths of propagation. A RAKE combiner 108 linearly combines these estimates to produce 
a higher confidence combined estimate of the original encoded data. Finally a decoder 110 extracts the original data 
transmission. 

15 [0009] The signal at the input of a channel estimation and compensation processor 106 is a modulated baseband 
signal distorted by propagation though the channel and by random noise. Data modulation may be based on any known 
methodology, but for purposes of the illustrative embodiment of the invention is taken to be quadrature phase shift key- 
ing (QPSK), which represents a data stream as a sequence of four-state symbols. 

[0010] Figure 2 illustrates the four QPSK symbol states. An in-phase baseband carrier signal I 200 is set to a + 1 
so state or a -1 state, the latter being 180 degrees away from a + 1 state. Similarly, a quadrature baseband carrier signal 
Q 202 is set to a + 1 state or a - 1 state, the latter being 180 degrees away from a + 1 state. I and Q signals are com- 
bined to produce a two-dimensional signal with the four states 204, 206, 208, and 210 shown in Figure 2. 
[001 1 ] For mathematical convenience, a time sequence of two-dimensional signals Eke those of a QPSK modula- 
tion may be represented as a complex function of time. Let complex time function z(t) be the originally transmitted 
25 QPSK-modulated baseband signal. Then the signal at the input of a channel estimation and compensation processor 
106 is: 

u(t) = c[t)z{t) + n(t) [Equation 1] 

30 where c(t) is the complex, time-varying channel gain and n{t) is random noise. If one can compute a good estimate of 
c(f) (denoted hereafter as c (/)) and n{t) is negligibly small or reduced to insignificance by averaging, then a good esti- 
mate of z(t) is: 



35 




[Equation 2] 



[001 2] Equation 2 can be rearranged to solve for c (f) and accordingly the approximate value of the channel distor- 
tion at any instant may be calculated from the values of the received signal and the original signal at that instant If 
40 known pilot symbols with values p are inserted in the original signal at specific points in time, then the original signal is 
known at those points and the value of the channel distortion at those points may be calculated from the value of the 
received signal at those points and the original pilot signal values. The "channel estimation" is then derived as an inter- 
polation of these sample values of the channel function at the pilot signal points. 

[001 3] Figure 3 shows a data stream 302 divided into ordered, consecutive blocks of M symbols each. The period 
45 of transmission of each block is T B . At the beginning of each block, one or more known pilot symbols are inserted in the 
data stream 302. 

[001 4] Figure 4 is a diagram of a channel estimation and compensation processor that uses inserted pilot symbols 
according to the method of Sampei and Sunaga (S. Sampei and T. Sunaga, "Rayleigh Fading Compensation lor QAM 
in Land Mobile Radio Communications'', IEEE Transactions on Vehicular Technology, Vol. 42, Mo. 2, May 1993, pp. 

so 137-147). An input analog baseband signal u{t) 400 is provided to and processed by a symbol synchronization module 
408 and a block synchronization module 406. Outputs of the symbol and block synchronization modules are processed 
to produce a pilot symbol sample clock 404 that controls an analog-to-digital sample circuit 402 (depicted with a switch- 
like symbol), which operates to cause digital values of pilot symbol signals to be sent to the complex multiplier 41 1 . The 
output of the complex multiplier 41 1 represents an estimate of the channel gain at the time of a particular pilot symbol. 

55 This estimate is sent to delay line 41 2. 

[001 5] The symbol synchronization module 408 also controls an analog-to-digital sample circuit 41 0 (depicted with 
a switch-like symbol) that operates to cause digital values of the input baseband signal u(t) 400 at each symbol time to 
be stored in a delay line 426. 
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[001 6] Considering Figure 3 and Figure 4 together now, a method for estimation of channel gain for the symbols in 
an arbitrarily chosen block 0 in the data stream 302 of Figure 3 will be descrtoed. With reference to the time frame 300 
in Figure 3, a symbol sample time 

'Km-WB + (53 7 fl [Equation 3] 



is the sample time of the mth symbol in the ton block To compute channel gain estimates for symbols in block 0, pilot 
10 symbols are sampled at t 10 . h o, and t 10 > ~ there is one pilot symbol at the beginning of each block, and the sym- 
bols from blocks -1 , 0, and 1 are used for channel estimation. For each trrpui pilot symbol sample u{t^ 0 ), complex mul- 
tiplier 41 1 computes 



c(t k .o ) = «('* * ) * Pk = c ('k o ) + '<<m o )/Pk (Equation 4} 



where p k is the known value of the pilot symbol at the beginning of data block k. These channel gain estimates are 
stored in the delay line 412. 

20 [0017] Complex multipliers 414, 416. and 418 are used to multiply the channel gain estimates in the delay line 412 
by appropriate weighting coefficients ct ,(m), a 0 (m), and a , (m), and the multiplier results are summed in complex adder 
420. [Values of the coefficients a are determined empirically according to known methods.] The resulting sum 422 (at 
output of adder 420) is therefore represented algebraically as : 



) = ^'o.o + « L ) = c(/ rt + — T 3 ) = or., (m)c{t_ l0 ) + or 0 (m)c(/ 0 u ) + a, (m M / ) 



30 [Equation 5] 



with T s defined to be the symbol period. 

Equation 5 constitutes a second-order Gaussian interpolation of the value of the estimated channel gain at the symbol 
35 sample time of the mth symbol in an arbitrarily selected block 0 based on the estimated values of the channel gain at 
the pilot symbol sample times stored in delay line 41 2. 

[0018] The inverse of resulting sum 422 is then computed by complex multiplicative inverter 424. A baseband sym- 
bol sample 428 is generated by delay line 426. where the delay is established to synchronize that baseband symbol 
sample with the computation of the estimated channel gain at the same sample time. Complex multiplier 430 then mul- 
40 tiplies the output of inverter 424 with symbol sample 428 to generate an estimate of an original data symbol 432 in 
accordance with Equation 2. 

[0019] Because the method of channel estimation just descrtoed inserts only one pilot symbol per data block and 

uses the pilot symbols from only three data blocks, the method is highly vulnerable to corruption by noise. 

[0020] In Figure 5, a channel estimation and compensation processor is depicted that uses inserted pilot symbols 

45 according to the method of Andoh, Sawahashi. and Adachi (See H. Andoh, M Sawahashi, and F. Adachi, "Channel 
Estimation Using Time Multiplexed Pilot Symbols for Coherent Rake Combining for DS-CDMAS Mobile Radio". (IEEE, 
PIMRC97, Helsinki. Finland, September 1-4, 1997), that method being Known as weighted multi-slot averaging 
(WMSA). With this method, an input analog baseband signal t/(f) 500 is processed by a symbol synchronization module 
408 and a block synchronization module 406. Outputs of the symbol and block synchronization modules are processed 

so to produce a pilot symbol sample clock 504 that controls an analog-to-digital sample circuit 402 (depicted with a switch- 
like symbol) to cause digital values of pilot symbol signals to be sent to complex multiplier 51 1 . The output of complex 
multiplier 511 represents an estimate of the channel gain at the time of a particular pilot symbol. This estimate is sent 
to delay line 512. 

[0021] The symbol synchronization module 408 also controls an analog-to-digital sample circuit 410 (depicted with 
55 a switch-like symbol) that operates to cause digital values of the input baseband signal u(t) 500 at each pilot symbol 
time to be stored in a delay line 526. 

[0022] With reference back to Figure 3, the estimation of channel gain for the symbols in an arbitrarily chosen block 
0 in the data stream 302 will be described for the WMSA method using the channel estimation and compensation proc- 
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essor of Figure 5. A symbol sample time, as depicted in Figure 3 and defined by Equation 3 above, is the sample time 
of the mth symbol in the kih block. Pilot symbols are sampled at f* jfn , where /c = -2. -1 , 0. 1 , 2. 3 and m = 0,1 , 2. 3 - 
i.e„ there are four pilot symbols at the beginning of each block, and the symbols from blocks -2, -1 . 0, 1 , 2, and 3 are 
used for channel estimation. For each input pilot symbol sample u[t k rr j, the complex multiplier 51 1 computes 

c(t km ) = u(t km )/p km = c(i km ) + n(t km )/p km [Equation6] 



io where k =-2, -1.0, 1, 2, 3, m = 0, 1. 2, 3, and p Km is the known value of a pilot symbol at position m in block A. These 
channel gain estimates are stored in the delay line 512. 

[0023] Complex averaging modules 534 are used to average the four channel-gain estimates associated with four 
. consecutive pilot symbols at the start of each block that are stored in the delay line 51 2. This averaging minimizes the 
error associated with the noise term of Equation 4. The noise-term error is, of "course, greater for the method described 
75 in conjunction with Figure 4. which does not use such averaging. 

[0024] Complex multipliers 514. 515, 516, 517, 518, and 519 are used to multiply the outputs of averaging modules 
534 by appropriate weighting coefficients a_2, <x. 1t clq, <x_i, a 2 , and a$, and the multiplier results are summed in complex 
adder 520. The resulting sum 522 (at the output of adder 520) is therefore represented algebraically as: 

^Co.o ) = *('«.- ) = Z -rZ £ <'* - > (Equation 7J 



25 

From this equation, it can be seen that the channel gain estimate c (t 0 Q ) is used for all symbol samples times in block 
0. regardless of the value of m. Also note that the a coefficients are independent of m. For the WMSA method, these 
coefficients are derived empirically by adjusting them during experiments to achieve the best performance for a partic- 
ular set of actual channel conditions. However, these coefficients are held constant during actual operation - i.e., they 

30 do not adjust dynamically to changing conditions 

[0025] The inverse of resulting sum 522 is then computed by complex multiplicative inverter 424. A baseband sym- 
bol sample 528 is generated by delay line 526, where the delay is established to synchronize that baseband symbol 
sample with the computation of the estimated channel gain at the same sample time. Complex multiplier 430 then mul- 
tiplies the output of inverter 424 with symbol sample 528 to generate an estimate of an original data symbol 532 in 

35 accordance with Equation 2. 

[0026] In summary, this WMSA method of estimating channel gain has better noise tolerance than the first method 
because of averaging of consecutive pilot symbols, but it uses static formulae for interpolation and does not dynamically 
adjust to changing channel conditions. 

[0027] An adaptation of the methodology for determining channel gain based on inserted pilot symbols, according 
<o to the method of the invention, is hereafter described in conjunction with the channel estimation and compensation 
processor depicted in Figure 6. An input baseband signal u(f) 500 is processed by a symbol synchronization module 
408 and a block synchronization module 406. Outputs of the symbol and block timing modules are processed to pro- 
duce a pilot symbol sample clock 504 that controls an analog-to-digital sample circuit 402 (depicted with a switch-like 
symbol) to cause digital values of pilot symbol signals to be sent to the complex multiplier 51 1 . The output of the com- 
45 plex multiplier 51 1 represents an estimate of the channel gain at the time of a particular channel symbol. This estimate 
is sent to the delay line 51 2. 

[0028] The symbol synchronization module 408 also oontrols an analog-to-digital sample circuit 41 0 (depicted with 
a switch-like symbol) that operates to cause digital values of the input baseband signal u(t) 500 at each pilot symbol 
time to be stored in a delay line 526. 

so [0029] The estimation of channel gain according to the method of the invention, for the symbols in an arbitrarily cho- 
sen block 0 in the data stream 302 of Figure 3. will now be described in relation to the functional elements of the channel 
estimation and compensation processor of Figure 6. A symbol sample time, as depicted in Figure 3 and defined by 
Equation 3 above, is the sample time of the mth symbol in the fcth block. Pilot symbols are sampled at t km where k »= - 
2. -1 . 0, 1 , 2, 3 and m = 0. 1 , 2, 3 - that is, there are four pilot symbols at the beginning of each block, and the symbols 

£5 from blocks -2. - 1 . 0, 1 . 2, and 3 are used for channel estimation.. For each input pilot symbol sample u(t k „i t channel 
gain estimates are computed according to Equation 6 by complex multiplier 51 1 . These channel gain estimates are then 
stored in the delay line 512. 

[0030] Complex averaging modules 534 are used to average the four channel gain estimates associated with the 
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four consecutive pilot symbols at the start of each block that are stored in the delay line 512. This averaging operates 
to minimize the error associated with the noise term of Equation 4. It should be understood, as well, that the use of more 
than four symbols per block would effect an even better cancellation of the noise error. Accordingly, while the preferred 
embodiment of the inventive methodology is based on the use of four symbols per block, that method also contemplates 
5 the use of a greater number of symbols. 

[0031] Complex multipliers 614, 61 5. 616. 617. 618, and 619 are used to multiply the outputs of averaging modules 
534 by appropriate weighting coefficients o_ 2 (m), cL^m). a 0 (m), a. 7 (/n), ct2(m), and c^m), and the multiplier results 
are summed in complex adder 520. The resulting sum 622 (at the output of adder 520) is therefore represented alge- 
braically as: 

10 

c(i 0 , m )^cit 0 o^mrj^i(i no ^^r a )^ Z^^XHi^ m ) [Equations} 

15 



The a coefficients are determined as follows 

(with q = a : 

20 



a(m) = -4K« 5 -yx«-3) [Equalion 9] 



a ( m ««-our-4x<r-3) [E ion l0 , 

30 |X 24 



35 



40 



a Q im) = -W 2 -1K9^)(9'3) [£quation , q 



a ,im) = 9L3±mg^xa& [Epuation 12 , 



a 2 (m )n -q(? 2 - 1 M^)<?- 3 ) [Equation 13] 

a 3 (m) = 9t3 ^3 [Equation 14] 

so [0032] Equations 8 through 14 constitute a fifth-order (six point) Lagrange polynomial interpolation of the value of 
the estimated channel gain at the symbol sample time of any symbol in block 0 from the estimated values of the channel 
gain at the pilot symbol sample times stored in delay line 512. According to the method of the invention, such interpo- 
lation may be made with a Lagrange interpolation of any order (r - 1), using the corresponding r blocks centered on 
block 0. and using any number P < N of pilot symbols at the beginning of each block. The resulting channel gain esti- 

55 mates are: 
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c« 0 .J = c(/„o +mr, ) = c(/. vil +£r 8 ) = Yr^-Jic(t im ) 



[Equation 15] 



where -0.5(/-2) <; k <> 0.5 r for even r. and -0.5<r-1) <; k < 0.5 (r-1) tor odd f. 
[0033] In this general case, the a coefficients are determined as follows 



10 



<withg = £) 



For even r: 
15 [0034] 



Mr 



20 



(0.5r - 1 + *)! (0.5r - *)'(</ - *) iJ 



fI(<7 + 0.5/--<j) [Equation lb} 



25 For odd r 
100351 

a * (lw) 3 (O.5r-0.5^)l[0.S(r-1)-/ f ll( Q ^) IT » * 0.5(M)-a] [Equation 17] 



(0036) With reference again to the channel estimation and compensation processor of Figure 6. the inverse of 
as resulting sum 622 is then computed by complex multiplicative inverter 424. A baseband symbol sample 528 is gener- 
ated by delay fine 526, where the delay is established to synchronize that baseband symbol sample with the computa- 
tion of the estimated channel gain at the same sample time. Complex multiplier 430 then multiplies the output of inverter 
424 with symbol sample 528 to generate an estimate of an original data symbol 632 in accordance with Equation 2. 

40 Conclusion 

[0037] A channel estimation methodology has been described that operates on one or more consecutive pilot sym- 
bols pet data block and incorporates Lagrange interpolation of the estimated channel gain between pilot symbols. The 
parameters of this interpolation adjust dynamically to changing channel conditions. 
45 [0038] Although the methodology of the invention, and illustrative applications of that methodology, have been 
descrfced in detail, it should be understood that various changes, alterations, and substitutions can be made therein 
without departing from the spirit and scope of the invention as defined by the appended claims. 

Claims 

so 

1 . A method for recovering a data signal having a plurality of data blocks wherein each said data block includes a plu- 
rality of coded data bits and a plurality of known pilot bits, said method comprising the steps of: 

estimating an amount of channel distortion in a received data signal based on comparison of received data sig- 
£5 nai relative to said known pilot bits; 

applying a Lagrange Polynomial Interpolation to said estimated channel distortion to provide a weighted esti- 
mate of channel distortion; 
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applying said weighted estimate of channel distortion to each of said data blocks to compensate for said chan- 
nel distortion. 

2. The method according to Claim 1, wherein said step of estimating channel distortion includes; 

5 

estimating said channel distortion from amplitude or phase information of said plurality of received pilot bits. 

3. The method according to Claim 2 wherein: 

io said estimate of channel distortion is made by computing an average amplitude value and an average phase 

value of said plurality of received pilot bits in said date blocks. 

4. The method according to Claim 1 wherein: 

15 said estimate of channel distortion is made by computing an average amplitude of said plurality of received 

pilot bits in said data blocks. 

5. The method according to Claim 2 wherein: 

20 said estimate of channel distortion is made by computing an average phase of said plurality of received pilot 

bits in said data blocks. 

6. The method according to Claim 1 wherein said step of applying a Lagrange Polynomial Interpolation includes the 
substep of: 

25 

applying a Lagrange Polynomial Interpolation of order R-1 to said estimate of channel distortion of R data 
blocks. 

7. The method according to Claim 6 wherein: 

30 

the weighting coefficients of said Lagrange Polynomial Interpolation are expressed by Equation 16 and Equa- 
tion 1 7 of the specification. 

8. In a receiving system for receiving a data stream comprised of a plurality of data blocks, said data blocks containing 
35 a plurality of known pilot bits and a plurality of coded data bits, a channel estimation and compensation processor 

comprising: 

distortion compensation means operative to determine a channel distortion parameter utilizing a Lagrange Pol- 
ynomial Interpolation weight factor algorithm: and 

40 

means for determining an estimate of channel distortion from said channel distortion parameter and applying 
said channel distortion estimate to said plurality of data blocks. 

9. The channel estimation and compensation processor of claim 8 further comprising: 

45 

a storage means for storing a portion of signal parameters of said plurality of pilot bits, said signal parameters 
being an estimate of channel distortion in related data blocks, and for storing at least one of said plurality of 
data blocks delayed in time by a predetermined time increment; 

so means to normalize each said estimate of channel distortion stored in said storage means, and means for 

applying said Lagrange Polynomial Interpolation computed weighting factors to each said normalized estimate 
of channel distortion; 

means for summing each said weighted normalized estimate of channel distortion, said sum being a estimate 
55 of channel distortion; and 

means for applying said estimate of channel distortion to said stored, delayed data block. 
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10. The channel estimation and compensation processor of Claim 9 wherein: 

said signal parameters are determined as an algebraic function of measured amplitude and phase values of 
said plurality of received pilot bits. 

1 1 . The channel estimation and compensation processor of Claim 9 wherein: 

said signal parameters are determined as an algebraic function of measured amplitude values of said plurality 
of received pilot bits. 

1 2. The channel estimation and compensation processor of Claim 9 wherein: 

_ said signal parameters are determined as an algebraic function of measured phase values of said plurality of 
pilot bits. 

13. The channel estimation and compensation processor of Claim 9 wherein: 



said signal parameters stored in said storage means correspond to a predetermined number, R, of said data 
blocks. 

20 

14. The channel estimation and compensation processor of Claim 1 3 wherein: 

said predetermined delay time increment is a first value if R is even and a second value if R is odd. 

25 1 5. The channel estimation and compensation processor of Claim 1 4 wherein: 

said first value of said predetermined delay time increment is equal to the integer value of (R/2+1) and said 
second value of said predetermined delay time increment is equal to the integer value of R/2. 

30 1 6. The channel estimation and compensation processor of Claim 13 wherein: 

said Lagrange Polynomial Interpolation weight factors are computed using a Lagrange Polynomial Interpola- 
tion of order (R-1) 

35 17. Trie channel estimation and compensation processor of Claim 1 4 wherein: 

said weight factors are expressed in Equations 16 and 17 of the specification. 

40 
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